Transitions to Line-Defect Turbulence in Complex Oscillatory Media 
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The transition from complex-periodic to chaotic behavior is investigated in oscillatory media 
supporting spiral waves. We find turbulent regimes characterized by the spontaneous nucleation, 
proliferation and erratic motion of synchronization defect lines which separate domains of different 
oscillation phases. The line defect dynamics is controlled by the competition between diffusion, 
which reduces line length and curvature, and phase-gradient-induced growth. The onset of each 
type of defect- line turbulence is identified with a non-equilibrium phase transition characterized by 
non-trivial critical exponents. 
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Two-dimensional reactive media with oscillatory dy- 
namics support a variety of spatio-temporal patterns in- 
cluding spiral waves. In the vicinity of the Hopf bifurca- 
tion, spiral waves are described by the complex Ginzburg- 
Landau equation (CGLE) [0J|]. Spiral waves can also 
exist if the local dynamics is complex-periodic or even 
chaotic While the basic features of such regimes 

are akin to those of the CGLE, a complete description 
of complex oscillatory media cannot be given in terms of 
the CGLE. For example, these media may undergo bifur- 
cations where the period of the orbit doubles at almost 
every point in space |4|J^. The rotational symmetry of 
spiral waves is then broken by the presence of synchro- 
nization defect lines where the phase of the local orbit 
changes by multiples of 27r. These defect lines have been 
observed in a super-excitable system and in experi- 
ments on the Belousov-Zhabotinsky reaction 0]. 

In this Letter, we study the fate of the synchroniza- 
tion defect lines as the system parameters are tuned to 
approach the domain where spiral waves have chaotic lo- 
cal dynamics. We show the existence of a new type of 
spatiotemporal chaos where the global temporal period- 
icity of the medium is broken by the spontaneous nucle- 
ation, proliferation and erratic motion of the defect lines 
separating domains of different oscillation phases. We 
describe the basic mechanisms governing the dynamics 
of the defect lines and provide evidence that the onset of 
each type of defect-line turbulence is a non-equilibrium 
phase transition with non-trivial critical exponents. We 
also study inhomogeneous media without spirals where 
line motion has a different nature and different scaling 
laws due to the absence of overall phase gradients. 

We study reaction-diffusion (RD) systems where the 
local kinetics is described by R(c(r, t)), a vector of non- 
linear functions of the local concentrations c(r, t). For 
simplicity, we assume that all species have the same diffu- 
sion coefficient D. While our considerations should apply 
to any RD system exhibiting a period doubling cascade 
to chaos, the calculations described here ||] were carried 
out on the Rossler model where Rx = —Cy — c^^ — 
Cx-\-Acy^ Rz = CxCz — Ccz-\-B. We investigate the behav- 
ior of the system as C increases, with other parameters 



fixed at A = 5 = 0.2 and D = 0.4. 




FIG. 1. The evolution of the medium in the r^2-line turbu- 
lent regime at C = 4.38. Square system of linear size L = 256. 
Panel (a) shows the Cz{r) field. The dashed line indicates 
the Qi defect Une. Panels (b), (c) and (d) are equally sepa- 
rated in time by 48 spiral revolutions and show the ^2(r) field 
color-coded by grey shades. The fine dashed lines in panel (b) 
mark the centers of the shock lines. 

Beyond the Hopf bifurcation, the system supports a 
spiral solution and infinitely- many spatially-blocked con- 
figurations of spirals coexist with spatially inhomoge- 
neous states without spirals. This multist ability is pre- 
served away from the Hopf bifurcation, even for C values 
corresponding to chaotic regimes. The blocked configu- 
rations form irregular cellular structures, similar to those 
observed in the CGLE. Cells are centered on spiral cores 
and their polygonal boundaries are delimited by shock 
lines where spiral waves from two neighboring cells col- 
lide. Fig. [^(a) shows a snapshot of the Cz{t) field for a 
simple configuration with two spirals. 

As C increases beyond the Hopf point, the Rossler 
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ODE system exhibits a period-doubling route to chaos 
followed by band-chaotic regimes intertwined with win- 
dows of periodic behavior. In spatially distributed media 
supporting spiral waves two period-doubling bifurcations 
take place at C 3.03 and C 4.075. These values are 
larger than the corresponding values for the ODE, 2.83 
and 3.86, respectively. These shifts in the bifurcation di- 
agram arise from the concentration gradients created by 
the spiral waves and their values depend on the spiral 
wavelength. In spatially inhomogeneous media without 
spiral waves the spatial gradients are small and the shifts 
of the bifurcation points are not detectable. The period 
doublings in media with spiral waves are necessarily ac- 
companied by the appearance of synchronization defect 
lines {ft curves) whose existence is required to reconcile 
the doubling of the oscillation period and the period of 
rotation of the spiral wave 

In the period-2 regime (3.03 < C < 4.075), a single 
type of synchronization defect line exists. These Q curves 
are defined as the loci of those points in the medium 
where the two loops of the period-2 orbit exchange their 
positions in local phase space and the dynamics is ef- 
fectively period- 1. The period-2 oscillations on opposite 
sides of the Q curve are shifted relative to each other by 
27r (a half of the full period). A medium with period- 
4 dynamics may support two types of synchronization 
defect lines corresponding to the two different possible 
types of loop exchange for a period-4 orbit. Across fti 
curves, which inherit properties of the Q curves, the two 
period-2 bands of the period-4 orbits exchange, leading 
to a ±27r phase shift across them. The 1^2 curves corre- 
spond to the exchange of loops within the two bands, a 
finer rearrangement of the local cycle. Along them, the 
dynamics is effectively period-2 and there is a An phase 
shift as the curves are crossed. 

Synchronization defect lines can be conveniently lo- 
cated by constructing scalar fields encoding the distance 
between loops of the period-4 orbit in local phase space. 
To this aim, we chose to take advantage of the regular 
succession of peaks in the local time series of c^, whose 
heights are in one-to-one correspondence with the vari- 
ous loops of the orbit. Calculating, at each point r, four 
such consecutive concentration maxima Ai (r) and order- 
ing them so that Ai{t) < A2{t) < As{t) < A4(r), one 
may construct the scalar fields ^i(r) = A4{r) — Ai{r) 
and ^2(1*) = ^4(1*) — ^3(r). In the period-4 case, ^i(r) 
and ^2(1*) take on fixed non-zero values at points in the 
medium away from spiral cores and shock lines and van- 
ish at points where the loop exchanges occur [0. Indeed, 
^i(r) decreases to zero on the Qi curves while ^2(1*) van- 
ishes on both the Qi and Q2 curves. In the following, 
we study the fate of the ft lines and use ^i(r) and ^2(1*) 
both to determine their length and to visualize them. 
The ^2(1*) field corresponding to Fig. 0(a) is shown in 
Fig. |l|(b). The thick vertical line connecting the spiral 
cores is an Qi curve, while the thinner lines are r^2 curves. 

On the shock lines, where the phase gradient vanishes, 
the local dynamics is approximately that of the Rossler 



ODE, and is thus always more advanced along the bifur- 
cation diagram. In particular, chaos first appears on the 
shock lines (for C 4.20). For C = 4.3, where most of 
the medium is still in the period-4 regime, two-banded 
chaos is seen on the shock lines (Fig. @(a)). These local- 
ized chaotic regions give rise to fiuctuations which may 
result in the creation of "bubbles" - domains delineated 
by circular ^2 curves (Fig. |l|(b)). For C < Cq^ 0^ 4.306, 
the bubbles are formed with a size smaller than a certain 
critical value and collapse shortly after their birth. As C 
increases beyond Cq^, the bubble nuclei begin to prolif- 
erate, forming large domains whose growth is limited by 
collisions with spiral cores or other domains. 

A typical life-cycle of a domain is illustrated in 
Fig. 0(b-d). The shock lines are nucleation sites of fto 
domains. Consider the two bubble-shaped nuclei indi- 
cated by arrows in panel (b) which were born in close 
proximity. In panel (c), they have coalesced, forming 
one rapidly-growing domain which then collides with its 
neighbor, leaving a shrinking internal domain (panel (d)). 
The contact of two Q2 lines always leads to their recon- 
nection and a reduction of their total length. The contact 
of Qi and Q2 lines leaves one Qi line. This event is ac- 
companied by a change of sign of the phase shift across 
the ft I line. 
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FIG. 2. Time series of Cz concentration maxima: a) 
C = 4.30, b) C = 4.42 and c) C = 4.7. Left panels show 
the local dynamics at a point in the cells, while the right pan- 
els show the dynamics on the shock lines. Time is in units 
of thousands of spiral revolutions. In panels (b) and (c) the 
crossings of the Cz maxima reflect the passages through the 
observation point of Q2 or Qi curves, respectively. 

The evolution of the size and shape of the Q2 line en- 
circling a domain is controlled by the balance of two com- 
peting factors: propagation along phase gradient directed 
toward spiral cores which results in line growth, and the 
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tendency of diffusion to eliminate curvature and reduce 
the length of defect lines. To investigate the interplay 
of these two factors, a series of simulations was carried 
out on a system without spiral waves, but with constant 
concentrations corresponding to those in the spiral core 
imposed along one pair of parallel boundaries. This ef- 
fectively creates "spiral core" boundaries emitting trains 
of plane waves which collide in the center of the system 
to form a straight shock line. In this case, the nucle- 
ated domains have very simple finger-like shapes, normal 
to the shock line, and consist of two straight segments 
with approximately semi-circular caps (Fig. ||). The do- 
main growth velocity normal to the core boundaries, v^, 
varies with the radius R of the arc of the growing tip 
diS v± = Vp — A/R where Vp :^ 0.126 is the velocity of 
a straight Q2 line parallel to the core boundaries, and 
A 0.658. The linear dependence of v± on 1/R shown 
in Fig. ^ confirms the effect of mean curvature on the 
velocity of Q line propagation. Since the width of small 
domains is approximately equal to 2i?, one can estimate 
the critical size that must be exceeded for domain pro- 
liferation. One finds Rc — 5.228, in good agreement 
with direct measurements from the observation of do- 
mains whose shape does not change with time. 
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FIG. 3. Dependence of the domain propagation velocity, 
v±, versus the curvature, 1/R, of the tip for C = 4.30. A 
typical finger-shaped Q2 domain is shown in the right panel 
along with the inscribed circle used to determine R. 

The transition to r^2-line-defect turbulence in media 
with spiral waves, which occurs around C = 4.3, changes 
the character of the local dynamics observed in the bulk 
of the medium. As the Q2 lines propagate, the associated 
loop exchanges result in an effective band-merging in the 
orbits of local trajectories so that they take the form of 
two-banded chaotic trajectories (Fig. ||(b)). Although 
the local trajectories retain their period-4 structure be- 
tween two passages of Q2 lines, the long-time trajectory 
cannot be distinguished from that of two-banded chaos. 
Thus, the global transition of the medium to defect-line 
turbulence can be characterized locally as intermittent 
band- merging. 

As the parameter increases further (C > 4.44) the lo- 
cal dynamics undergoes prominent changes. It fails to 
exhibit a period-4 pattern in the intervals separating line 
defect passages, and consists instead of four-banded or- 
bits whose bands grow in width with increasing C and 



merge at C ^ 4.7. Together with this permanent band- 
merging, spontaneous nucleation of Q2 bubbles occurs in 
the bulk. These chaotic Q2 lines are the loci of medium 
points where the two chaotic bands of the local orbit 
shrink and a thick "period-2" orbit is formed. As C in- 
creases, their width decreases, and for C ^ 4.8 the Q2 
lines cease to exist as well defined objects. 

While the local dynamics changes continuously to four- 
and subsequently two-banded chaos, another transition, 
mediated by moving Qi lines, takes place. At C = 4.557, 
the shock regions, where the local dynamics exhibits one- 
banded chaos, begin to spontaneously nucleate bubbles 
delineated by Qi lines. As C increases, the newly-born 
domains begin to proliferate. The qualitative features of 
this transition are similar to those of the r^2-line turbu- 
lence transition: the dynamics of the lines encircling 
domains is controlled by the factors discussed above and, 
considering shape of the long-time local phase space tra- 
jectories, it can be associated with intermittent band- 
merging (Fig. §|(c)), leading to one-banded local chaotic 
orbits. 

As the parameter C increases even further, beyond 
C > 5.0, the local trajectories in the bulk of the medium 
exhibit complete band- merging to one-banded chaos. In 
this regime, no defect lines can be identified and further 
increase in C does not result in any qualitative changes. 
However, spiral waves continue to exist, signalling the 
robustness of phase synchronization in this amplitude- 
turbulent regime. [O 
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FIG. 4. Critical properties of the two defect-line mediated 
phase transitions observed, (a) Time series of p2 at C = 4.34 
for the configuration shown in Fig. ^ (L = 256, two spirals; 
time is in units of thousands of spiral revolutions), (b) Vari- 
ation of p2 with C [o - L — 256, two spirals; • - L — 256, 
four spirals; x - L = 512, two spirals), (c) pi versus C for 
an inhomogeneous system without spiral waves, L — 256. (d) 
Same as (b) but for the Qi-line transition. The solid lines are 
power- law fits. 
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We now focus on the two onsets of synchronization 
defect hne turbulence. The line density, pi{t) = 
ii{t)/y/S, where S is the surface area of the medium 
and £i{t) the instantaneous total length of Qi lines, can 
serve as an order parameter to characterize these tran- 
sitions. Above each transition threshold, and as long 
as the corresponding defect lines continue to exist, the 
balance between line growth and destruction results in 
a statistically stationary average density p^, while pi{t) 
fluctuates. Thus the time series of p2{t) above the first 
threshold shown in Fig. ^(a) exhibits high-frequency, low- 
amplitude fluctuations attributed to the birth and death 
of nuclei in the shock regions, as well as large-amplitude 
oscillations with long correlation time. This suggests 
that the proliferation of domains and their destruction 
through coalescence occurs cooperatively. This is con- 
firmed by the fact that, for both transitions, the order 
parameter goes continuously to zero as C decreases to- 
ward the threshold. (In Fig. ^(d) the pi density does 
not vanish below threshold because the contribution from 
the stationary line shown in Fig. |l| has not been re- 
moved.) The data fall on curves with power-law forms, 
~pi{C) ^ {C — CqJ^% the signature of continuous phase 
transitions. The critical values are found to be Cq^ 
4.306 and Cq^ 4.557, while the critical exponents are 
— 0.25 and /3i 0.49. Finite-size effects usually ac- 
company critical point phenomena as correlation lengths 
diverge near threshold. Here, the finite-size to consider 
is the typical size of the cells composing the spiral wave 
structure. Strictly speaking, pi is not an intensive quan- 
tity because line-defect motion is constrained to occur 
between the network of shocks (where they nucleate) and 
the spiral cores, and this area varies from one spiral con- 
figuration to another. However, the data in Fig. ||(b,d) 
show that pi{C) depends weakly on the cell size. 

A number of conclusions can be drawn from the sim- 
ulation results. The two transitions exhibit significantly 
different scaling properties, a remarkable fact given that 
the mechanisms at play appear to be the same in both 
cases. This difference may arise from the fact that the 
Qi line transition takes places when Q2 lines still exist in 
the medium. 

The zero-spiral-density limit is singular since the tran- 
sitions observed in the medium without spiral waves are 
different from those described above. In this case, the 
onset of defect-line nucleation occurs at the same critical 
values Cq.. However, in the absence of large-scale phase 
gradients, the entire medium behaves like the shock re- 
gions separating spiral wave cells, defect lines do not grow 
and the increase of 'pliC) arises essentially from the en- 
hanced nucleation rate. This leads to a different form of 
the onset of line turbulence (cf. Fig. ||(c) for the behav- 
ior of p2{C)) characterized by different critical exponents 
{PI = 1.22, (3^ = 0.53). These values are difficult to esti- 
mate because of fluctuations in the fields not associated 
with fully developed Q lines. The fact that and are 



significantly different from /3i and (32 supports the con- 
clusion that the character of the transitions is different 
in the spiral and spiral-free systems. 

To our knowledge, there is no equilibrium equivalent 
of these phase transitions, nor were their non-equilibrium 
analogs reported previously. The line-defect phase tran- 
sitions may constitute a special class of non-equilibrium 
critical phenomena since, in this form of spatiotemporal 
chaos, it is the dynamics of one-dimensional synchroniza- 
tion defects that breaks the global temporal periodicity 
of the medium. 

Finally, the observations of synchronization defect 
lines in a number of excitable systems |^,|^ demonstrate 
their independence of any particular reaction mecha- 
nism. Therefore, the phenomena presented in this Letter 
should be observable not only in chemical media exhibit- 
ing period-doubling but also in a much broader class of 
systems. For example, they may exist in the cardiac mus- 
cle where complex-excitable dynamics and spiral waves, 
necessary prerequisites for the emergence of synchroniza- 
tion defects, have been established experimentally 
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